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Abstract 


Can a violent process like sudden reheating after supercooling at the onset 
of a first-order QCD transition improve the possibility of primordial black 
hole formation? Underdensities reheat earlier than overdensities, there is a 
short period of huge pressure differences, hence ffuid acceleration. Density 
perturbations on scales far below the Hubble radius A <C Rn get an ampli¬ 
fication which grows quadratically in wavenumber, the amplifications at the 
horizon scale are small. Primordial black hole formation cannot be sufficiently 
amplified by the QCD transition unless the initial spectrum is fine tuned. 
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I. INTRODUCTION AND CONCLUSION 


The transition from a quark-gluon plasma to a hadron gas in the early universe took 
place at a temperature T* ~ 150 MeV. Recent results of lattice QCD indicate a hrst-order 
phase transition for the physical values of the u, d, s quark masses [|I|. The mass inside 
the Hubble horizon at the QCD transition is ~ IMq. Crawford and Schramm first 
proposed that black hole formation at the QCD transition could account for dark matter 
today. Recently Jedamzik [Q proposed to identify such primordial black holes (PBH) with 
the macho’s (massive compact halo objects) observed by microlensing [Q. He pointed 
out that the formation of PBH’s should be particularly efficient during the QCD transition 
due to a signihcant decrease in the effective sound speed during the transition. Schmid, 
Schwarz and Widerin showed that during a hrst-order QCD transition the sound velocity 
Cs = {dp/dpYj"^ vanishes (for wavelengths much larger than the bubble separation), there are 
no pressure gradients, no restoring forces, and preexisting cosmological density perturbations 
go into free fall. They computed the amplihcation of linear cosmological perturbations due 
to the QCD transition. The amplihcation is very large for scales A far below the horizon, 
amplihcation oc k'^ with n = 1 in the bag model, n = 3/4 in a ht to lattice QCD data. There 
is no amplihcation for superhorizon modes in accordance with a general theorem. Near the 
horizon the hrst peak (above the incoming spectrum of cosmological perturbations) has a 
height (amplihcation factor) of only 1.5 for the ht to the lattice QCD data and 2.0 for the 
bag model. The position of the hrst peak is at Mrad ~ O.lhM©, which corresponds to 
A/i?H ~ 0.6 at the end of the hrst order QCD transition. It was concluded that primordial 
black hole formation is unlikely, because the amplihcations of perturbations due to the QCD 
transition are large only far below the Hubble scale. 

In this paper we ask, whether a violent process like sudden reheating after supercooling at 
the beginning of the QCD transition could enhance the amplihcation of preexisting density 
perturbations and improve the possibility of black hole formation at the QCD transition. 

The amount of supercooling (T* — Tsc)/T* = r/ is proportional to the dimensionless ratio 
(cr3/V(/Ty^), where a = surface tension, I = latent heat, see e.g. 0]. Data are available 
only for quenched lattice QCD (gluons only, no quarks) l/T^ ~ 1.40(9) and a very small 
surface tension u/T^ ~ 0.0155(16). This gives a very small supercooling of ?7 ~ 10“^ [0. In 
view of the uncertainties in the lattice QCD determinations and the absence of data for I and 
a from QCD including three physical quarks we consider the possibility of a larger surface 
tension and/or a smaller latent heat giving more supercooling. For illustrative purposes in 
our hgures we used rj = 1/10. 

The amplihcation of perturbations due to sudden reheating at the beginning of a hrst 
order QCD transition occurs because reheating happens at a certain supercooling temper¬ 
ature Tgc. Hence underdensities reheat earlier than overdensities. (For subhorizon scales 
Newtonian time is applicable.) During this short period of time, there are spatial variations 
of pressure of 0{ri) leading to huge pressure gradients (compared to the the small preexisting 
cosmological pressure gradients) and huge huid accelerations, hence preexisting overdensi¬ 
ties get a very violent but very short compressional impulse. The huid velocity ehectively 
jumps proportional to the wave number k. For scales far below the Hubble radius this jump 
becomes the dominant contribution to the huid velocity throughout the transition. In the 
bag model IQ, which we shall consider here for simplicity, the resulting amplihcation of 
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density perturbations for scales far below the horizon grows as i.e. with one additional 
enhancement factor oc k due to sudden reheating. For horizon and superhorizon scales 
general relativity must be used. 

The general relativistic time-evolution requires the choice of a foliation of space-time, 
and for linear perturbations this dehnes a gauge. The hypersurface of sudden reheating is 
a slice of constant energy density, and we use the corresponding gauge, the uniform density 

(UD) gauge, to evolve through sudden reheating. In this gauge the relevant variables are 
continuous at sudden reheating. To evolve before and after reheating the UD gauge is not 
suitable, because it is singular in the subhorizon limit, instead we use the uniform expansion 

(UE) gauge, which is nonsingular, and which is defined by requiring that the fundamental 
observers, who are at rest on a slice, have a uniform Hubble expansion. In the UE gauge 
the relevant variables are discontinuous at sudden reheating, and these discontinuities are 
obtained in a very simple way by using the gauge transformation between the UD and the 
UE gauges. In the superhorizon limit the state variables in UE gauge do not jump. Moreover 
we prove the theorem that phase transitions do not affect at all the laws of evolution for 
6 = 6p/p of the growing mode (i.e. growing in relative importance) in UE gauge in the limit 


(<5ue) 


growing mode 


= const 
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This translates into the ’conservation law’ (‘^ue) growing mode ~ const, where 4 A (puE = 
by using the energy constraint. Because of this theorem the superhorizon 
modes during the QCD transition remain unaffected in the sense of Eq. (|ID. Near the 
horizon the first amplification peak gives a factor 2.2 for a supercooling of 10% in the bag 
model, the additional enhancement due to sudden reheating at the horizon scale is about 
10%. 


The enhancement of black hole formation at the QCD transition is discussed in the last 
section. For standard models of structure formation without a large tilt, the amplitudes 
are not big enough to produce a cosmologically relevant amount of black holes 0. A 
tilted spectrum could be hne tuned to produce black holes at the QCD scale. But this 
tilted spectrum would need a break just below the QCD scale in order not to overproduce 
smaller black holes, a second fine-tuning. We conclude that the QCD transition enhances 
the probability of black hole formation, but for an observable amount of black holes today 
the preexisting spectrum would have to be fine tuned around the QCD scale, and the major 
effect would not be due to the QCD transition. 


II. SUPERCOOLING AND SUDDEN REHEATING 

The quark-gluon plasma, photons and leptons are tightly coupled via strong, electromag¬ 
netic and weak interactions, V/H S> 1, and make up a single perfect (i.e. dissipationless) 
radiation fluid at scales A > see [0. The pressure for perfect fluids with negligi¬ 

ble chemical potential only depends on temperature. In the bag model p], which we will 
consider for simplicity, the quark-gluon plasma (QGP) is described by Pqqp = Pqgp ~ 
i.e. 


3 




(2) 


Pqgp - ~ 

where g* is the effective number of relativistic helicity states, and B is the bag constant. 
Above the phase transition temperature we include u, d quarks and gluons in the quark- 
gluon plasma and 7 , e, /r, p's in the photon-lepton fluid. Below the phase transition we model 
the radiation fluid as a hadron gas (HG) of massless pions, tightly coupled to the photons 
and leptons. 

In the cosmological expansion the temperature in the bag model drops oc 1/a, see Fig. |I|. 
To derive this result we use s = dp/dT, which comes from the Maxwell relation for the 
free energy, hence Sqgp = Sqct> oc T^. We further use that the entropy in a comoving 
volume is conserved, s oc 1/a^, and finally obtain T oc 1/a. The energy density is given by 
p = Ts — p from the first law and homogeneity, hence p^Qp = Pqqp + B. It follows that 
the sound velocity (Cg)QGP = {dp/dT)^^^ = 1/3 and that p/p < 1/3. The bag constant is 
determined by the transition temperature T* and /S.g* = 5'qgp~5'hg viapQGp(T*) = phg(^*), 
B = /90)Ag*T^. Similarly the latent heat is given via I = T^^As, I = (27r^/45)A5f*T/. 

For supercooling, T < T*, the free energy density / = —p(T) for the QGP is higher 
than for the HG, and therefore the QGP is not the equilibrium state. However the fluid 
supercools in the metastable quark-gluon phase until the formation and growth of hadronic 
bubbles starts reheating the universe at Tgc < T*. We give a short overview of bubble 
nucleation and sudden reheating, a more detailed discussion of bubble nucleation can be 
found e.g. in 0. Assuming homogeneous bubble nucleation (no ’dirt’) the probability to 
form a bubble is proportional to exp(AS'), where AS is the change in entropy by creating 
a bubble. AS is determined by the work required to form a (spherical) bubble, —T*AS' = 
AF = (pqgp — Phg)^-R^ + crdvri?^, where a is the surface tension of a hadronic bubble. 
Only for large enough bubbles {R > Rent) free energy is gained. The probability to form a 
hadronic bubble with critical radius per unit volume and unit time is given by 

/(r)=/„exp(-^^) , (3) 

with AFcrit = 167r(T^/[3(pHG ~Pqgp)^]- For small supercooling, p = 1 —T/T* -C 1, we obtain 
(Pqgp ~ Phg) ~ (T — T^)A{dp/dT) = —pT^As = —pi, where I is the latent heat. Therefore 
the probability to form a critical bubble per unit volume and unit time can be written as 

I^T^e^p{-A/p^), (4) 

with A = 167 r(j^/( 3 /^T*). A numerical prefactor to T/ would be irrelevant for our purposes. 

The surface tension a and the latent heat I are the crucial parameters for Tgc. Data are 
available only for quenched lattice QGD (gluons only, no quarks) 0: l/T^ ^ 1.40(9) and 
a very small surface tension cr/T/ 0.0155(16). There are no values for unquenched QGD 
available yet. Using the results from quenched lattice QGD we find A = 2.9 x 10“®. 

The bubbles grow most probably by weak deflagration EHIl. and the latent heat re¬ 
leased from the bubbles is distributed into the surrounding QGP by acoustic waves and by 
neutrinos, whose mean free path is 10“^i?H, of the order of the bubble separation. This re¬ 
heats the QGP to T* and prohibits further bubble formation. The supercooling temperature 
fraction p can be estimated by the schematic case of one single bubble nucleated per Hubble 
volume and per Hubble time. 
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This rough estimate gives a hnal supercooling fraction which is only 20% too low compared 
to the more realistic case of one bubble nucleated per cm^ per 10“® of a Hubble time, see 
0. The time needed for the supercooling is given by Atgc/tn = = 0(10“^). In the 

figures we took rj = 10“^ for illustrative purposes. At the maximal supercooling the universe 
suddenly reheats up to the QCD transition temperature T*. The sudden reheating implies 
a jump in pressure [p] = p* — Psc- 


\P] 

P + P 
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up to linear order in p. This jump in pressure depends on the maximal supercooling and 
can be a huge effect. In the case of the cosmological QCD transition it is ~ 100 times larger 
than preexisting CODE normalized perturbations. On the other hand, the jump in entropy 
density is small since it has to be quadratic in p due to the second law, [s] /s = |p^. 

The reheating is treated as a sudden process happening at the maximal supercooling, 
see Fig. |l|. This is justihed since the reheating time is ~ Atgc, which follows from 

the following consideration. Bubbles present at a given time have been nucleated typically 
during the preceding time interval Atnuci = I/{dI/dt). Using the relation between time and 
supercooling rj, drj/dt = Sc^/tn, we hnd 


AWiAh = vViGAcl) = 0(10-®). 


(7) 


The corresponding bubble nucleation distance is a few cm ||l4|,[15 


During the reversible part of the first order phase transition thermal equilibrium between 
the quark-gluon plasma and the hadron gas is maintained. We showed in [§ that pressure 
gradients and the isentropic sound speed (for wavelengths A much larger than the bubble 
separation), c* = {dp /must be zero during a hrst-order phase transition of a fluid with 
negligible chemical potential (i.e. no relevant conserved quantum number). The sound speed 
must be zero, because for such a fluid the pressure can only depend on the temperature, 
p{T), and because the transition temperature T* has a given value, it cannot depend on 
any parameter, hence p(T*) = p* is a given constant, and Cg = 0. During the entire QCD 
transition the sound speed stays zero and suddenly rises back to the radiation value c* = 
1/v^ after the transition is completed. In contrast pressure varies continuously and goes 
below the ideal radiation fluid value p = p/3, but stays positive. 


III. AMPLIFICATION OF SUBHORIZON PERTURBATIONS 

The evolution of linear cosmological perturbations is hrst analyzed on subhorizon scales, 
where Newtonian concepts for space and time (for an expanding radiation fluid) are appli¬ 
cable. The perturbations on subhorizon scales get an amplihcation factor which is quadratic 
in wave number k for H. 

The maximal supercooling is reached when bubble nucleation becomes efficient and the 
radiation fluid suddenly reheats, i.e. it happens on a hypersurface of constant temperature. 
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T = Tsc, hence uniform energy density. On this hypersurface of locally sudden reheating, 
Srh, pressure jumps uniformly, see Fig. 

At a given Newtonian time, overdensities have higher temperature than underdensities 
and therefore reheat later. This time delay (lapse of time) of the actual reheating of a fluid 
element on Srh from the average Newtonian time of reheating, is denoted by At(x) and 
shown in Fig. |]for an overdensity of one perturbation mode k. The time delay At(x) follows 
from the condition of uniform energy density on the surface of reheating, p(x, f + At(x)) = 
uniform. 


ex, 


i(-) 


Af(x) 


dp 

dt 


(-) 
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( 8 ) 


e(x) = 6p denotes the energy density perturbation at a given time, (—) means hxed time 
immediately before the short period during which the various fluid elements reheat. Inserting 
the continuity equation for the FRW background dp/dt = —3i7(p+p) at the time of maximal 
supercooling gives 


Af(x) 
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?,H p + p 


(9) 


One can apply the same argument for the perturbations and the background just after 
reheating. The time delay of Srh is the same whether evaluated before or after the hyper¬ 
surface, and this condition gives the discontinuity equation for e. 


e 
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p + p 

“ p + p 

p + p 


0 . ( 10 ) 


p is continuous at the actual hypersurface Srh, but the effects of the fluid expansion is 
different in the two media. This generates an effective discontinuity in e. 

We now integrate the evolution equations for the relativistic fluid over the very short 
time period of order At when reheating happens at some x. We split p and p into a 
homogeneous term and an inhomogeneous term, p(x,f) = p{t) -|- e(x, t), p(x, t) = p{t) + 
7 r(x, f). In distinction to usual perturbation theory the pressure inhomogeneities compared 
to the background energy density p are of order AT/T*. The pressure inhomogeneities are 
the dominant terms, the driving terms, all remaining terms can be neglected during the 
short time interval when reheating occurs at various x. In this approximation the continuity 
equation (energy conservation) and the Euler equation (momentum equation) read. 


dte{'x,t) = —3if7r(x, t) (Ha) 

cIiS(x,t) = -V7r(x,f), (11b) 


where S is the momentum density and S = (p -|- p)v for a radiation fluid. 

The pressure inhomogeneity 7r(x, t) is a sequence of step functions with step size [p], 
while V7r(x, t) is a sequence of Dirac delta functions. It is necessary to distinguish the 
unperturbed pressure p{t), which makes a jump [p] at t from the average pressure 

p{t), which does not make a jump at t If is very convenient to dehne the pressure 
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perturbations as vr(x, t) = p(x, t) — p(t), therefore vr(x, t) 7 ^ p(x, t) — p{t). The dominant 
pressure inhomogeneity is 


7 r(x, t) = \p] {9(t - tun - At(x)) - 9(t - tun)}, (12) 

with [p] = — p^~\ Integrating the continuity equation Eq. (|lla|) over the short period 

of huge pressure differences gives the jump condition for e, 


[e(x)] = [p] 3HAt{x) 


[P] 


6(X) 


p + p 


(-) 


(13) 


in agreement with Eq. 0 ). 

Integrating the Euler equation Eq. ( |1 1 b] ) over the short period of huge pressure differ¬ 
ences, gives the jump condition for the momentum density S, 


[S(x)] 


[p] / dtS{t 


r 1 1 „ e(x) 
\p] —V- 

p + p 


tRH - At(x))V(Af) 
(-) 


(14a) 

(14b) 


The discontinuity in S is a local discontinuity at the actual hypersurface of reheating (pro¬ 
duced by the 6 functions of Eq. ( |llb| ) and Eq. (|14a|) ). It is a sudden impulse imparted on 
each fluid element. 

In the evolution equations Eqs.( |lla|) , ( |llb| ) different Fourier modes are coupled during 
the short period of reheating. Therefore the standard procedure to evolve each k mode 
separately is not applicable at reheating. However the resulting jump conditions obtained 
by working in x space, Eqs. (p!3D , (|14b| ), are valid for each k mode separately. 

On subhorizon scales the discontinuity in e is independent of k (for a preexisting Harrison- 
Zel’dovich spectrum), the discontinuity in S is proportional to fcphys and is the dominant 
effect in the subhorizon limit. 

The evolution of density perturbations before and after sudden reheating can be com¬ 
puted in a Jeans analysis generalized to a relativistic fluid for each Fourier mode separately, 
see p. The preexisting subhorizon perturbations before sudden reheating are acoustic oscil¬ 
lations, see Fig. H. In the bag model the sound speed Cg = 1/a/ 3 and the amplitudes for 
e and |v^S| are equal until sudden reheating. At ^rh, e and S jump according to Eqs. (0), 
(PEI) . After sudden reheating, during the reversible part of the QCD transition, the sound 
speed Cg = 0 , i.e. the restoring force in the acoustic oscillations vanishes as hrst shown in 
and density perturbations go into free fall. Since the QCD transition lasts less than 
a Hubble time, gravity is negligible for subhorizon scales and density perturbations have 
approximately constant velocity, e{t) = — (t — %h)VS*'^^ This gives an amplihcation 

which is linear in fcphys and proportional to After the QCD transition is completed, 

one has acoustic oscillations again and the total amplihcation factor for k ^ H is 


A 


out 


k , k 


= 7-|7-cos(pi„ 
Ain ki /C 2 


sm(pin| 


(15) 


</ 9 in is the phase of the incoming acoustic oscillation at ^rh- The scale ki = v^/Attrans ~ 
depends on the duration of the reversible QCD transition, Attruns- This contribution is due 
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to the free fall and starts just below the Hubble scale, see Fig. |^. k 2 depends on the amount 
of supercooling, ^2 = ifRH2A/3/?7. The amplihcation factor is quadratic in k for (X/Rh) < rf. 

Jump conditions like Eqs.(|I^, (|14b|) must be used whenever there is a jump in the 
background quantities (as opposed to a jump only in as in the case without supercooling 
and sudden reheating). Only the background pressure is allowed to jump not the energy 
density, according to the continuity equation. Due to the second law of thermodynamics a 
fluid can only heat instantaneously, sudden cooling is impossible. 


IV. GENERAL RELATIVISTIC ANALYSIS 


A slicing of space-time (with space-like hypersurfaces on which time is dehned to be 
constant) is needed to formulate time-evolution equations in general relativity. We restrict 
ourselves to time-orthogonal foliations. For linear cosmological perturbations a choice of a 
foliation already hxes a gauge. For sudden reheating there is one special slice, the hypersur¬ 
face of reheating Srh, which is a surface of constant energy density. It is a hxed-time slice 


in uniform density (UD) gauge, see |T6 


We use the uniform density gauge, the gauge adapted to our problem, to evolve through 
sudden reheating. We now show that in UD gauge the relevant variables are continuous (do 
not jump). The hypersurface of sudden reheating, Srh, has dehnite extrinsic curvature A, 




and intrinsic curvature at each point on Srh, hence these quantities are continuous at 

reheating. The momentum density S in UD gauge also stays continuous at reheating due to 
the momentum constraint of general relativity, DjK\ — DiKR = SnCSi, see e.g. El- The 
energy density is uniform in space and continuous in time in UD gauge. 

The UD gauge is not suitable for evolving before and after reheating, because the UD 
gauge is a singular gauge in the subhorizon limit, iJ/fcphys —^ 0, see below. Therefore we shall 
solve the dynamics before and after sudden reheating using the uniform expansion gauge 
(UE), which is nonsingular both in the superhorizon and subhorizon limits and in which 
the physical dynamics will turn out to be most transparent (in contrast to the UD gauge). 
However in the UE gauge the relevant variables are discontinuous (jump) at sudden reheat¬ 
ing. The uniform (Hubble) expansion gauge is dehned by requiring that the fundamental 
observers, who are at rest on the slice, M(obs) = n(S), have uniform Hubble expansion, i.e. 
the perturbation of the mean extrinsic curvature of S, 6 triF* (S) = k, is zero, kue = 0. 

The geometrical properties of a hypersurface S in the longitudinal sector are given by 
the following variables, Bardeen 1989 The perturbation of the trace of the extrinsic 
curvature 6 triFj(E) = k, the traceless part of the extrinsic curvature (shear aj of normals), 
which is determined by x through o'j (E) = — (d''dj — y, and *^^^i?(S) = —4 A ip. The 

lapse of time between the perturbed hypersurfaces is (1 -|- a). The connection between the 
geometric variables and the gauge potentials (in the longitudinal sector with time-orthogonal 
coordinates), 


= — (1 -|- 2a) dt^ -|- a{tY [Jp (l -l- 2ip) -|- 2didj'y\ dx^dx\ (16) 

is as follows: y = a^y, k, = — Ha) — Ay. 

The state of the fluid (before or after sudden reheating) on a hypersurface of hxed time (in 
some gauge) is given by the following variables: The energy density perturbation Sp = e and 













the momentum density S, which can be written as S = 'W'lp in the longitudinal sector (sector 
of scalar perturbations). p{p) is given by the bag model before and after the irreversible 
reheating. There are no anisotropic stresses for our perfect fluid (QCD, 7 , leptons). The 
total set of variables in the longitudinal sector is (k, y, (p), a. 

For the dynamical equations in UE gauge we choose to take the evolution equations of 
the fluid, = 0 , which are the continuity equation and (in the longitudinal sector) the 

3-divergence of the Euler equation, 

dte = —3H {e + Ti) — A'lp — 3H(p + p)a (17a) 

dtip =— Ti — {p-\-p)a^ (17b) 

with TT = 5p. The system of dynamical equations is closed by Einstein’s i^Qg-equation, 

(A + 3if)a = 47 rG(e + 37 r), (17c) 


together with the equation of state. This set of equations has exactly the same structure as 
the Jeans equations: Two hrst-order time-evolution equations for the fluid state variables 
(e,'^)uE, whose initial values can be chosen free of constraints, supplemented by an elliptic 
equation for ctue, which plays the role of the gravitational potential. The lapse a is the 
only geometrical variable needed to solve the dynamical evolution in UE gauge. The other 
geometric perturbations, i.e. the shear y and the intrinsic curvature 93 , do not appear in our 
closed set of evolution equations, Eqs. ( tl7a| - |17c|) . If desired, ip and y can be computed from 
the state variables via the energy and momentum constraints, /\p = —dvrGe, Ay = —\2ttG'iP 
on any UE slice. 

To make the gauge transformation from the uniform expansion gauge (independent vari¬ 
ables e, V’) to the uniform density gauge (independent variables -ip, k) we need the gauge 
transformation formulae for (e, " 0 , k) 


e — e =—3H{p + p)At (18a) 

ijj —'ll) =—{p + p)At (18b) 

k- K = -{3H - kl^y^)At , (18c) 


where the variables with tilde correspond to the new gauge and At is the lapse of time from 
the old to the new hypersurface, see Fig. ^ 

We now derive the discontinuity conditions at sudden reheating for the uniform expansion 
gauge. We use the state variables (e, V’)ue before reheating and use Eq. ( |18a|) with 
euD = 0 to obtain 


Af(-) 


^UE 

3H{p + pG)y 


(19) 


where (—) means fuE immediately before reheating. The momentum potential -^ud and the 
extrinsic curvature kud follow from Eqs.( |18b|) , (|18c|) . 


^UD 


i’isE 


^UE 

3H 


k\jd — -{3H^ ^ - /Cpi^yJ 


^UE 

3H{p + p(-)y 


( 20 a) 

( 20 b) 
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'j/'UD and KuD stay continuous at sudden reheating. We make the analogous gauge trans¬ 
formation (from UE to UD) after sudden reheating. From the kud equation we get the 
discontinuity condition for e, 


3H - fcgi,,. 

P+P 


Cue 


= 0 , 


( 21 ) 


where the jump in the background quantity H = —47rG(p -|- p) is the given input, and the 
jump in e is the output. In the subhorizon limit [^/{p + p)] = 0 and i.e. 


we 


recover the results of the analysis using Newtonian geometry, Eq. (|T^. In the superhorizon 
limit e does not jump, [e] = 0, but ^ This means that the coordinate regions 

%E < and tuE > ^RH cover certain regions of space-time doubly and other regions not 
at all. Applying the b'uD equation before and after reheating we obtain 


[V'ue] - ^ [^ue] ■ 


( 22 ) 


The resulting modihcation of the perturbation amplitude due to the QCD phase transition 
(transfer function) is shown in Fig. 

On superhorizon scales (e,'ib)uE do not jump. An even stronger statement can be made 
for 5 ue = ^ue/p: The evolution of the growing mode (i.e. growing in relative importance) 
in (5 ue is given by 


(■luEjgrowing = const 



(23) 


for scales |i^|} and p/p = w < 1. Written in this form the growing mode of the 

state variable ^ue is manifestly unaffected during phase transitions or any other changes in 
the equation of state. The proof goes as follows. We take the ^phys l-^l limit of the R^q- 
equation, —?,{p + p)a = (e-t-Svr), and we insert this into the continuity equation, Eq. (|17aj) , 
and into the Euler equation, Eq. (|17b|) , 

dte + 2He = (24a) 

dti) + SH'ijj = -e. (24b) 

3 

We note that the pressure gradients have dropped out. For the mode which will turn out to 

be the growing mode the two terms on the left-hand side of the Euler equation do not cancel, 

and e is relevant in the Euler equation, hence xjj = O {e/H). Therefore the term is negligible 

in the continuity equation for the growing mode if k^^^^ 77^, c?t(lne) = —2c}t(lna). Hence 

e = const and using the Friedmann equation 5 ue = ^ue/p = const (/Cphys/77)^. For 

the decaying mode the two terms on the left-hand side of the continuity equation do not 

cancel, and xjj is relevant in the continuity equation, hence e = O (^physb’/^)- Therefore 

the e term is negligible in the Euler equation for the decaying mode, ^^(In'^) = —3c?t(lna). 

2 

Hence -0 = const a~^, and the continuity equation gives e cx a~^ f a~^dt. With a cx it 

follows that this is indeed the decaying mode for w = p/p < 1. The law, Eq. (|23|), can be 
translated (using the energy constraint) into the ’conservation law’ (^Pue) growing mode ~ const. 
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This is consistent with the ’conservation law’ [|T^ Cgrowing mode = const for /Cp^ys -C where 


Deruelle and Mukhanov |T^ have analyzed jump conditions (matching conditions) in 
zero shear (ZS) gauge, also called longitudinal gauge, where the shear of the normals to 
the equal time hypersurfaces is zero by dehnition. They focussed on superhorizon physics. 
The matching conditions are much more complicated in zero shear gauge than in uniform 
expansion gauge; no superhorizon discontinuities in UE state variables, complicated discon¬ 
tinuity conditions in ZS in Ref. 0. This is connected to the fact that the zero shear gauge 
is singular in the superhorizon limit (Bardeen 1988 jl^), and one particular aspect of this 


is that Atzs.uo/'^^UE.UD oc (Tf/fcphys) —>■ oo in the superhorizon limit. 

Finally we show that the uniform density gauge is singular in the subhorizon limit. For 
a Harrison-Zel’dovich spectrum the subhorizon density contrast f^UE = e/p = 0(10 is 
independent of fc, the state variables (e, '0 )ue agree with (e, -0) in a Jeans analysis. From 
Eq. (|T^ and the subhorizon limit of Eq. (pOb|) we obtain k\]J)/H = 0(/Cpj^yg(5uE/-ff^) ^ oo. 


Eq. 


gives Suo/euE = 0(kphys/iJ) ^ oo. 


V. BLACK HOLE FORMATION AT THE QCD TRANSITION ? 


Black holes form in a radiation dominated universe if the density contrast of a top hat 
perturbation inside the Hubble radius is in the range 1/3 <Jh< 1[^|- For an observable 
amount of IMq black holes today, i.e. = ^(1); ^Fe fraction (5 of energy density 

converted to black holes at the QCD transition must be Cl(aQCD/aequaiity) ~ 10“®. For a 
gaussian distribution of density fluctuations the fraction (3 of pen at the time of formation 
is given by. 




Pbh 


Ptot v^27rj ] 


1 


1 /*^ * ' 2^2 
rms ^rms 


(26) 


Without any enhancement from the QCD transition this requires hrms ~ 0.06 H- The 
sudden reheating at the onset of a hrst-order QCD transition leads to huge amplihcations of 
density perturbations on scales far below the Hubble horizon. But at the horizon scale the 
QCD transition gives enhancement factors of 2.0 and 2.2 for the bag model without resp. 
with sudden reheating for a supercooling of 10%, Fig. an additional enhancement of 10%. 
For lattice QCD without sudden reheating the enhancement factor is 1.5 ^ in our linear 
perturbation treatment. This indicates a corresponding reduction in the required preexisting 
perturbation spectrum at the solar mass scale. Cardall and Fuller used a qualitative 
argument of Carr and Hawking and the bag model and also obtained a factor 2 reduction 
in the required preexisting perturbation spectrum. These QCD factors of 1.5 or ~ 2 are so 
modest that a preexisting Harrison-Zel’dovich spectrum with COBE normalization is very 
far from giving a cosmologically relevant amount of black holes [jl^ • One would have to put 
in a hne-tuned tilt (n — 1) ~ 0.36 to get the desired amount of black holes. However, this 
tilted spectrum would overproduce primordial black holes on scales which are only a factor 
50 below the Hubble radius at the transition. Therefore a break in the preexisting spectrum 
below the QCD scale would be required, a second hne tuning. 
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We conclude that the QCD transition with or without sudden reheating enhances the 
probability of black hole formation, but the preexisting spectrum needs to be hne tuned 
around the QCD scale, and the major effect would not be due to the QCD transition. 
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FIGURES 



FIG. 1. The evolution of temperature T as a function of the scale factor a in the bag model. 
For illustrative purposes the supercooling is AT/T* = 10“^. 



FIG. 2. A space-time diagram at the time of sudden reheating: The supercooled (2) and the 
reheated (4) epoch are separated by the hypersurface of reheating Srh (3). On Srh temperature, 
hence pressure, jumps uniformly. At(x) denotes the lapse of time compared to the average time, 
tRH) for a certain fluid element to reach the hypersurface of reheating. For subhorizon wavelengths 
tRH is the average Newtonian time for reheating, for horizon and superhorizon wavelengths (general 
relativistic case), t = tRH is a hypersurface of constant time in a given gauge, e.g. in the uniform 
expansion gauge. 



14 





FIG. 3. The evolution in conformal time rj of the density contrast (Jrad = ^/p) and the 
velocity (V’rad = • S/{kp) = (p + p)v/p) for the radiation fluid in uniform expansion (Hubble) 

gauge. At sudden reheating, the fluid velocity jumps. During the reversible part of the QCD 
transition, marked by the 2 vertical lines, the velocity stays approximately constant and the density 
contrast grows linearly. The incoming amplitude is normalized to 1. 
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FIG. 4. The modifications of the radiation fluid amplitude Aout/^in (transfer function) due 
to the QGD transition in the bag model with a supercooling of AT/T* = 10“^. On the horizontal 
axis the wavenumber k is represented by the RAD mass contained in a sphere of radius Ti/k. Mh 
is the mass inside the horizon at the QGD transition. 
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